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The Hydrodynamics of Superfluid Turbulence (HST) describes the flows (or 
counterflows) of Hell in the presence of a chaotic set of vortex filaments. 
The HST equations govern both a slow variation of the hydrodynamic vari- 
ables due to dissipation related to the vortex tangle and fast processes of 
the first and second sound propagation. This circumstance prevents effec- 
tive numerical simulations of the problems of unsteady heat transfer in HeU. 
By virtue of a pertinent multi-scale perturbation analysis we show how one 
can eliminate the fast processes to derive the evolution equation for the slow 
processes only. We then demonstrate that the long-term evolution of a tran- 
sient heat load of moderate intensity obeys the nonlinear heat conductivity 
equation. The second example of the methods developed is investigation of 
unsteady processes in the second sound resonator. The latter is frequently 
used for study of nonstationary behavior of vortex tangle, just by monitor- 
ing of the quality factor behavior. This procedure however is wrong when 
characteristic times of processes are comparable (or smaller) than the time 
constant of resonator. We show how to extract the correct information on 
the vortex line density (VLB) dynamics with use of procedure we developed. 

PACS numbers: 47.32.Cc, 47.37.-hq, 67.40.Vs., 05.10.Gg. 



1. INTRODUCTION AND SCIENTIFIC BACKGROUND 

It is appreciated that in most cases the flow (or counterflow) of Hell 
occurs in the presence of a chaotic set of vortex filaments, the so called 
superfluid turbulence. To describe various problems both for applications 
and for experiments it is necessary to use equations of HST. These equations 
unify classical equations of hydrodynamics of HeH and the famous Vinen 
equation describing the dynamics of the vortex tangle. The HST was widely 
used to describe various hydrodynamics as well as thermal phenomena in 
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Hell in presence of vortex tangle (see the review paper 1^. 

The set of HST equations is extremely cumbersome, therefore, there is 
no wonder that to achieve quantitative results one is bound to turn to nu- 
merical methods. Yet, numerical simulation of nonstationary flows of Hell 
faces one serious obstacle. The point here is that slow variation of hydrody- 
namic variables is accompanied by the fast processes related to propagation 
of the first and second sounds. If one is interested in the slow evolution, 
particular details of sound propagation become completely irrelevant, yet 
requiring rather extensive numerical resources. 

Thus, it seems attractive to try to get rid of the fast modes by a pure 
analytical procedure. In the present paper we realize effective separation 
of the slow from the fast modes using multi-time asymptotic perturbation 
techniques (See, for example, book by Nayfeh^^. That is quite universal 
powereful method, its application, however, depends on concrete statement 
of the problem. In the present work we consider two particular examples. 
In the first example we follow how the initially hyperbolic-type equations of 
HST (describing wave type of heat transfer) lead to purely parabolic-type 
nonlinear heat transfer equations. Another example of the method is the 
study of unsteady processes in the second sound resonator. 



2. THE NONLINEAR HEAT CONDUCTIVITY EQUATION 

Let us consider the problem of one-dimensional heat exchange in the 
tube filled by Hell. We assume that a heater is switched instantly on and 
the thermal front starts out propagating into undisturbed bulk of helium. 
Initially, of course, the propagation of the heat should be realized by the sec- 
ond sound mechanism. But gradually, due to quantum vorticity the second 
sound attenuates and degenerates. The mechanism of the heat exchange 
should change. In order to ascertain both this new mechanism and the 
correspondent law we apply multi-time asymptotic perturbation techniques. 

The set of equations of HST for dimensionless velocity of the normal 
component , dimensionless temperature T' and dimensionless square root 
of the vortex line density (VLD) G = yj C/ Coo ( for details and notations 
see the review article of Nemirovskii and Fiszdon ^ is reduced to the form: 

^ + ^ - shG'v ^ + ^ -0 m 

dt' ^ dx' - dt' + dx' - °' 

dG a ,a, Sh , 



dt' A{T)psPn 'P' 2 



(GV^ - G% (2) 
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with the following dimensionless variables: 
t = -t', X = Dx', Vn = VnoV^, T = ^T', C = CooGh (3) 

C2 aTC2 

It is easy to see that the dynamics of the hydrodynamic variables is 
specified by the dimensionless criterion the Strouhal number Sh = D/c2Td, 
(D-length of channel, Td decrement of attenuation of counterflow due to 
vortices). The Strouhal number is defined as the ratio of the counterflow 
decrement to inverse time that takes the heat pulse to cross the channel. In 
the case of heat load of moderate intensity, of order of a few W/cm?, the 
quantity Sh is of order of 10^ — 10^, i.e. much greater than unity. 

For times of order L/{c2Sh) the heat pulse propagates according to 
wave-like equation. Our goal now is to clarify what will happen at later 
times. For that purpose we have to eliminate the fast processes. According 
general ideas of multi-scale asymptotic perturbation theory (see e.g.Nayfehl^ 
we look for a solution to the set of the HST equations (EQ)-© in the form of an 
asymptotic series. Following this method, we introduce different time scales. 

tQ = t'; t'l = et'; t'2 = eH' , where e = l/Sh « 1. We 

look for a solution to the set of the HST equations ©-((H) in the form of an 
asymptotic series = V^{x', t'o, t[,t'2)+eVl{x' , t'o,t[,t'2)+e^V^ix' , t'o,t[,t'2) + 
.... and similarly to T' and G. The simple rule takes place: 

The next step in study of the slow evolution of the heat pulse consists in 
substituting the multi-time scale series into equations of HST (P)-©- Gath- 
ering terms of the same order of magnitude with respect to e we come up with 
a chain of equations leading to divergent (secular) solutions. Cancelling step 
by step these secularities we then obtain a hierarchy of equations of different 
orders in parameter e, governing different stages of the evolution. Omitting 
details of computations we concentrate on the first order in e. In the first 
order, equations for normal velocity and temperature variations read 

dV^ _ dT^ ^ -0 

dt[ dx'"^ ' dt' dx' 

Excluding quantity Vq from the set of equations written above we obtain 

dt'^ dx' \ dx' ) ' ^ ' 

Relation © coincides with the widely used nonlinear heat-conductivity 
equation derived by Dresnerl^by a pure phenomenological model. We stress, 
that unlike the Dresner derivation, our work shows that this relation is just 
the long-time asymptotic limit of the full set of equations of HST. 
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Fig. 1. (a) Excess attenuation of the resonance waves measured by author 
of 1^, (b) The same quantity computed with use of asymptotic perturbation 
techniques. We have to stress that this behavior is just the transient re- 
sponse of the resonator, and change of amplitude does not follow (directly) 
the evolution of VLD. 



3. MEASUREMENTS OF THE VORTEX LINE DENSITY IN 
THE SECOND SOUND RESONATOR 

The second problem which we study with use of multi-scale asymptotic 
perturbation techniques is scrutinizing a very popular experimental method: 
measuring vortex line density with a second sound resonator. Although that 
method is standard and used in many works, we, for the sake of definiteness, 
will consider the recent experiments made by Hilton and Van Sciverl^, They 
performed direct measurements of quantum turbulence induced by second 
sound shock pulses in a wide channel. In Fig. la. there is depicted the excess 
attenuation of second sound in the resonator placed across the channel, along 
which the intense thermal pulse is propagating. It was assumed that thermal 
pulse created the quantum turbulence which lead to the excess attenuation 
of second sound. Assuming further that the the excess attenuation at,(t) 
is connected with the amplitude of resonance wave A(t) by the well-known 
Hall-Vinen relation ay{t) = ay{0){A{0)/A(t) — 1) , and expressing a^(t) 
via the VLD L(t) the authors drew a conclusion about the VLD evolution, 
coinciding (up to factor) with the curve depicted in Fig. la. 

It is easy to see, however, that this evolution looks to be quite unrea- 
sonable. Indeed, according to current view on VLD dynamics, the quantity 
L{t) should grow only in presence of relative velocity. But the pulse dura- 
tion was 0.75 ms , whereas L(t) increases during 5-^7 ms. The origin of 
that discrepancy is obvious. The Hall-Vinen relation based on the quality 
factor is valid only for very slowly changing external conditions (duration of 
thermal pulse) in comparison with the time constant of the resonator. Here 
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the condition was inverted and the Hah-Vinen relation was inapphcable; i.e. 
the amphtude A{t) evolved independently from the quality factor evolution. 
In other words the curve depicted in Fig.l. demonstrates reaction of the 
resonator rather than the true behavior of VLD L(t). 

Thus the direct procedure of measuring of VLD based on the quality 
factor is not correct and the question appears whether it is possible at all to 
extract information about VLD L(t) evolution. The answer is positive and 
corresponding analysis is based on the general topic of our work i.e. on the 
fact that excess attenuation of the second sound due to quantum turbulence 
changes slower than duration of main thermal pulse. Then we are in position 
to apply the multi-scale asymptotic perturbation techniques. 

Let us describe very briefly the corresponding computations. Taking 
the spacial Fourier transform to Eqs. ^ we arrive at the following equation 
for variations of the temperature T' in the standing second sound wave 

^|:,,(0^+.V = .=T„e" (6) 

Here Tg is the pumping amplitude, uj is the resonance frequency of 
standing wave, p = a;, + a^(^) is the attenuation coefficient consisting of 
background part due to background vorticity and variable part ) due 
to passing thermal pulse. The latter depends on the slow time ^ = et, where 
e is small parameter, characterizing, say, attenuation coefficient/frequency. 
The corresponding analysis demonstrates that the solution is (as usual) a 
combination of general and particular solution. The general solution is just 
natural oscillations, attenuating proportional to 
exp[— (1/2) / p{(,)dS,]. As for particular solution, we seek it in a form: 



r;,,, = r,{t, 6)e-*, ^ = [G{C, e) - i^M + ^(C, e) (7) 



Substituting ((T)) into © and comparing coefficients near r\{t^ e)e*'^* and e"^* 
we obtain relations for G and H. Representing, further, quantities G and 
H in form of series in e, and equating coefficient near the same powers of 
e we obtain recurrent equations for successive determination of both G'(^) 
and H{S^). Restricting ourseves to the first order we obtain the following 
equation to determine the complex amplitude of particular solution, which 
is just time dependent amplitude in a resonator measured in experiment!^ 

-dt= 2 ^ + (' )^ 

It is seen that amplitude of second sound wave in resonator is equal to 
the pumping amplitude multiplied by the quality factor only in limit of very 
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long time. Generally that amplitude obeys a more complicated law. Let us 
discuss how to extract the information concerning true behavior of vortex line 
density. Here we have performed the following procedure. We supposed that 
the Vinen equation for evolution of VLD can be indeed extrapolated to very 
high value of applied counterflow velocity. Then we varied coefficients aand /3 
of that equation in order to select best values to fit the curve depicted in curve 
of Fig. la. In Fig. lb. is a curve which corresponds to values a ~ 0.13 ayjand 
/3 « 0.16 where ayjand /3yi are values offered by Vinen. Evolution of 
the VLD corresponding to these values is depicted in Fig. 2. 

Fig. 2. Plausible behavior of VLD re- 
sulting in response of resonator de- 
picted in Fig.L This curve was ob- 
tained as solution of Vinen equation 
with parameters obtained by fitting 
procedure. One can see, that VLD 
increases only during thermal pulse 
duration 0.75 ms. 

4. CONCLUSION 

By use of multi-scale perturbation analysis it has been shown how the 
wave-like behavior of the heat pulse described by initially hyperbolic equa- 
tions gives way to the parabolic-type non-linear heat conductivity equation. 
By the similar method we study unsteady processes in the second sound 
resonator. 
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